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lO . ABSTRACT 

J Context. Dust particles evolve in size and lattice structure in protoplanetary disks, due to coagulation, fragmentation and crystalliza- 

Ph , tion, and are radially and vertically mixed in disks due to turbulent diffusion and wind/radiation pressure forces. 

^0 ' Aims. This paper aims at determining the mineralogical composition and size distribution of the dust grains in planet forming regions 

of disks around a statistical sample of 58 T Tauri stars observed with Spitzer/IRS as part of the Cores to Disks (c2d) Legacy Program. 
Methods. We present a spectral decomposition model, named "B2C", that reproduces the IRS spectra over the full spectral range 
(5-35 /jm). The model assumes two dust populations: a warm component responsible for the 10/jm emission arising from the disk 
O ■ inner regions (< 1 AU) and a colder component responsible for the 20-30/jm emission, arising from more distant regions (< 10 AU). 

5-H ' The fitting strategy relies on a random exploration of parameter space coupled with a Bayesian inference method. 

j/j ' Results. We show evidence for a significant size distribution flattening in the atmospheres of disks compared to the typical MRN 

C^ , distribution, providing an explanation for the usual flat, boxy 10/jm feature profile generally observed in T Tauri star spectra. We 

reexamine the crystallinity paradox, observationally identified by Olofsson et al. (2009), and we find a simultaneous enrichment of 
the crystallinity in both the warm and cold regions, while grain sizes in both components are uncorrelated. We show that flat disks tend 
to have larger grains than flared disk. Finally our modeling results do not show evidence for any correlations between the crystallinity 
>^ , and either the star spectral type, or the X-ray luminosity (for a subset of the sample). 

\l ■ Conclusions. The size distribution flattening may suggests that grain coagulation is a slightly more effective process than fragmenta- 

\l ' tion (helped by turbulent diffusion) in disk atmospheres, and that this imbalance may last over most of the TTauri phase. This result 

>M ' may also point toward small grain depletion via strong stellar winds or radiation pressure in the upper layers of disk. The non neg- 

>—' I ligible cold crystallinity fractions suggests efficient radial mixing processes in order to distribute crystalline grains at large distances 

r~^ . from the central object, along with possible nebular shocks in outer regions of disks that can thermally anneal amorphous grains. 
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1 . Introduction in disks, and therefore constitute a very important ingredient in 

any planet formation theory. 



d ' The mid-infrared spectral regime probes the warm dust grains 

located in the planet forming region (typically 1-10 AU for a The dust grains that are originally incorporated into pro- 
classical T Tauri disk). At these wavelengths, the young disks toplanetary circumstellar disks are essentially of interstellar 
are optically opaque to the stellar light, and the thermal emis- nature. They are thought to be particles much smaller than 
sion arises from the disk surface, well above the disk midplane. 1 yu m, and mostly comp osed of silicates or organic refracto- 
Single-aperture imaging of disks in the mid-IR suffers from both ries. iKemper et alJ (l2005h placed an upper limit of 2.2% on the 
a relatively low spatial resolution compared to optical/near-IR amount of crystalline silicates in the interstellar medium (here- 
telescopes, and from poorly extended emission zones. On the after ISM), which suggests an amorphous lattice structure for 
other hand, silicates have features due to stretching and bend- the pristine silicates in forming protoplanetary disks. In the very 
ing resonance modes which make spectroscopy at mid-IR wave- early stages of the disk evolution, the tiny dust grains are so cou- 
lengths of very high interest and feasible with current instru- pled with the gas that grain-grain collisions occur at sufficiently 
mentation. Silicates are indeed the most abundant sort of solids low relative velocity to allow the grains to coagulate and grow. 

This results in fractal aggregates that will tend to settle toward 
the disk midplane as their mass increases. From simple theo- 
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forces on the grains in a laminar disk, one can show that /im- 
sized grains at 1 AU should settle in the midplane in less than 
10^ years in a classical T Tauri disk (Weidenschilling 1980). In 
fact, the grains are expected to settle even faster as their mass 
increases in the course of their journey to the disk midplane. 
Because the T Tauri disks are a few million years old, this would 
suggest that the inner disk upper layers should be devoided of 
//m-sized or larger grains in the absence of turbulence and grain 
fragmentation. This is a prediction that can be tested observa- 
tionnaly, especially in the infrared where spectroscopic signa- 
tures of silicate grains are present. 

The Spitzer Space Telescope, launched in August 2003, had 
a sensitivity that surpassed previous mid-IR space missions by 
orders of magnitudes until the cryogenic mission ended in May 
2009. As part of the "Core to Disks" (c2d) Legacy survey (Evans 
et al., 2003), more than a hundred T Tauri stars were spec- 
troscopically observed, to confirm or invalidate some of the 
predictions concerning dust processing and grain dynami cs in 
protoplanetary disks. In iKessler-Silacci et aP (l2006l l2007h . and 
lOlofsson et al.l (l2009l) . we showed that most of the objects dis- 
play silicate emission features arising from within lOAU from 
the star This allowed both a classical study of grain coagulation 
and a comprehensive statistical analysis of dust crystallization in 
planet forming regions of di sks around young solar analogs. 

In lOlofsson et al.l ( |2009|) , we showed that not only the warm 
amorphous silicates had grown, but so had the colder crystalline 
silicates as probed by their 23jum emission feature. In fact, the 
emission features in IRS spectra are very much dominated by 
micron-sized grains in the upper layers of disks, pointing toward 
vertical (turbulent?) mixing of the dust grains to compensate for 
gravitational settling, together with grain-grain destructive frag- 
mentation in the innermost regions of most protoplanetary disks 
to compensate for grain growth. This equilibrium seems to last 
over several millions of years and be to indep endent of the spe- 
cific star forming region dOliveira et al.ll2010l) . Winds and/or ra- 
diation pressure, in complement to these processes, can act to 
remove a fraction of the submicron-sized grains from disk atmo- 
spheres, and may thus contribute to the transport of crystalline 
si licate grains toward th e outermost disk regions (see discussion 
in lOlofsson et alj|2009l) . 

Crystalline silicates appear to be very frequent in disks 
around T Tauri disks and in regions much colder than their 
presumed formation regions, suggestin g efficient outward ra- 



dial transport mechanisms in disks ( Bouwman et al.l l2008t 
lOlofsson et al]|2009l: I Watson et alj|2009al) . Therefore.lhe deter- 
mination of the composition and size distribution of the dust 
grains in circumstellar disks is one of the keys in understand- 
ing the first steps of planet formation as it can trace the dust 
fluxes in planet-forming disks. Detailed modeling of the lOjum 
silicate e mission feature h a s alrea dy be en performed for HAEBE 
disks by iBouwman eTall ( 1200 ll) and Ivan Boekel et al.l ( l2005h . 
Mineralogical studies of dust in disks around ve r y low mass 
stars a n d brown dwarfs we re al so led by Apai et al. (l2005l) . lRiaa 
(120091) . iMerin et al.1 (l2007h and lBouvet all C2008.) . The two last 
studies introduced a novel method to fit the spectra over the en- 
tire IRS spectral range, which allows to decompose the spec- 
tra into two main contributions at different temperatures. This 
in turn allowed to compare the degrees of crystallinity and the 
grain sizes in two different disk regions. The latter composi- 
tional fitting approac h is further supported by the analysis by 
lOlofsson et al.l (l2009l) who showed that the energy and frequen- 
cies at which crystalline silicate features are seen at wavelengths 
larger than 20 //m are largely uncorrected to the amorphous 
IQ/im feature observational properties. 



We present in this paper an improved version of t he composi- 
tional fitting method used in lMerin et al.l (l2007h and lBouv et al.l 
(2008). We apply the mod el to a subsample (58 stars) of high 
SNR spectra presented in lOlofsson et alJ (120091) to derive the 
dust content in disks about young solar analogs. The method re- 
lies on the fact that the IRS spectral range is sufficiently broad 
so that the regions probed at around IQfixn and between 20 and 
30 yum do overlap only parti ally. This is illustrated for instance in 
iKessler-Silacci et al] (|2006|) who show that for a typical T Tauri 
star, a factor of 2 increase in wavelength (lOyum — > 20;um), 
translates into a factor of 10 (1 AU —> 10 AU) into the regions 
from which most of the observed emission arises. An additional 
goal of our model is therefore to search for differential effects 
in crystallinity and grain sizes between the warm and slightly 
cooler disks regions which may be indicative of radial and/or 
vertical dependent chemical composition and grain size, due for 
instance to differential grain evolution and/or grain transport. 

We develop in Sec.|2]our procedure to model IRS spectra of 
Class II objects, and present the tests for robustness of this pro- 
cedure in AppendixlAl The results for 58 objects are presented 
in Sec. [3] and they are discussed in terms of dust coagulation and 
dust crystallization in Sec.|4] Sec.|5]summarizes the implications 
of our results on disks dynamics at this stage of evolution and 
critically discuss the use of the shape and strength of the amor- 
phous lOyum silicate feature as a grain proxy. 

2. Spectral decomposition with the B2C method 

We elaborate in this section a modeling procedure whose goal is 
to reproduce the observed IRS spectra, from about 5 to 35 jum, 
in order to infer the composition and size of the emitting dust 
grains. This is achevied by using two dust grain populations at 
two different temperatures. We will refer to these two popula- 
tions as two different "components". The first component aims 
at reproducing the features at around 10 /im, its temperature is 
generally around T„ ~ 300 K (warm component hereafter). The 
second component aims at reproducing the residuals, over the 
full spectral range and its temperature is colder, T^ ~ lOOK 
(cold component hereafter). Each of these two components is 
described by several grain compositions, including amorphous 
and crystalline silicates, and sizes as detailed below. Bayesian 
inference is used to best fit the IRS spectra and to derive un- 
certainties on the parameters, hence the name of the procedure, 
"B2C", which stands for Bayesian inference with 2 Components. 

2. 1 . Theoretical opacities and grain sizes 

To reproduce the observed spectra, we consider five different 
dust species. The amorphous species include silicates of olivine 
stoichiometry (glas sy MgFeSi04, density of 3.71 g.cm"^, op- 
tical indexes from iDorschner et alj 1 1995b . silicates of pyrox- 
ene stoichiometry (gl assy MgFeSiOe, density of 3.2 g.cm"^, 
IDorschner et al.l Il995l). and silica (arnorpho us quartz, density 
of 3 33 g.cm~^ iHenning & Mutschkelll997h . For the Mg-rich 
(see lOlofsson et al.l l2009l) crystalline s pecies, we consider en- 



statite (MgSiOs, density of 2.8 g.cm ^ llaeger et al.lll998l) and 
forste rite (Mg2Si04, density of 2.6 g.cm"^, IServoin & Pirioul 

119731) . 

Following prev ious papers (e.g. lBouwman et ani200lLl2008l 
iJuhasz et al.l 120091) . the theoretical opacities ka of amorphous 
species are computed assuming homogeneous sphere (Mie the- 
ory), while the DHS theory (Distribution of Hollow Spheres, 
iMin et al1l2005h is employed for the crystalline silicates in or- 
der to simulate irregularly-shaped dust particles. To limit the 
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Fig. 1. Blowup of the opacities (in units of cm^.g ') used to model the IRS spectra. The solid line is for 0.1 yum grains, dashed line 
for l.Syum grains and dotted line for 6.0yum grains. 



number of free parameters in the model, we consider three spec- 
troscopically representative grain sizes radii, which are 0.1 //m, 
1.5 yKm and 6.0 yum. These three grain sizes are supposed to best 
mimic the behaviour of ver y small grains, intermediate- sized 
grains and large grains (e.g. iBouwman et ani200lL l2008l) . We 
limit ourselves to the two smallest sizes, 0.1 and l.Sjum, for 
the crystalline species. The first reason for this choice is that 
large crystalline grains present a high degeneracy with large 
amorphous grains, as 6.0/zm-sized grains are mostly feature- 
less. Therefore large crystals can be used by the procedure in- 
stead of large amorphous grains, thereby introducing a bias to- 
ward high crystallinity fractions. The second rea son is that, ac- 
cording to crystallisation models (e.g. lGailll2004l) . one does not 
expect to produce such large pure crystals via thermal anneal- 
ing. Grains more likely grow via collisional aggregation of both 
small crystalline and amorphous material (e.g. Min et al. 2008). 
Following these two considerations, we decided not to include 
6.0 Ai m-sized crysta l s. This ch oice i s in line with previous works 
from ISargent et al] (l2009l) or iRiaJ (l2009l) . The fifteen opacity 
curves (5 grain compositions, 3 grain sizes) used in this paper 
are displayed in Fig.[T] including 6.0jum-sized crystals to show 
the degeneracy with large amorphous grains. 



2.2. The B2C model 

The B2C model elaborated in this paper assumes that the IRS 
spectrum can be fitted by considering a continuum emission. 



and two main components, a warm and a cold one, essen- 
tially responsible for the IQjjm and 20-30 /im emissions, respec- 
tively. The two compo nent approach is supported by the work of 
lOlofsson et al.l (l2009h who show that disks usually have inho- 
mogeneous compositions with respect to the dust temperature 
(crystalline features being more frequently detected at long than 
at short wavelengths, the so-called crystallinity paradox), and 
that the emission features at around 10 //m are essentially uncor- 
rected with those appearing at wavelengths larger than 20yum. 
This approach is also supported by the T Tauri model described 
in Kessler-Silacci et al. ( 2006) where they show that the emis- 
sion at lOyum is arising from regions within ~ 1 AU while the 
flux at 20 /im is arising from within ~ 10 AU. Furthermore, we 
first used the simpliest solution, with only one component, with- 
out successful results. Any compositional method that aims at 
fitting IRS spectra over the full spectral range should therefore 
be able to handle inhomogeneous disk compositions. The two 
component approach developed in this paper is a simple attempt 
to go in that direction. This modeling strategy has already proven 
sufficient to provide adequate fits to IRS spectra, from 5 to 35 /im 
(Merin et al. 2007 ; Bouy et al. 2008 ). In this paper, we improve 
both the model an d th e fitting stra t egy or iginally developed in 
iMerin et al.l (l2007l) and lBouvefaD (l2008l) in order to apply the 
decompositional method to a large number of objects. 

The first step of the modeling approach is the estimate of the 
continuum to be subtracted to the observed IRS spectrum before 
performing a compositional fit. This could be done using a ra- 



Olofsson et al.: Spectral decomposition of C2D/IRS spectra of TTauri disks 



diative transfer code (e.g. iMerin et al.1 120071: iBouv et al.ll2008h . 
but given the number of objects (58) to be analyzed, and given 
the objectives of the paper which is oriented toward searching 
for trends thanks to the analysis of a large sample, obtaining a 
satisfying model for every object is not a manageable task. We 
instead adopt a continuum built by using a power-law (/l") plus 
a black-body at temperature Tcont, to make our problem more 
tractable: 



Sz96 



KpjA" + KbijByiTcont) 



(1) 



where Kpi and /Tbb are two positive constants. The power law 
index, a, is determined independently from the compositional 
procedure. It represents the mid-IR tail of the emission from both 
the star itself and from the inner rim of the disk. Its value is es- 
timated by fitting the IRS spectrum from its shortest wavelength 
(about 5fim) up to blue foot of the lOjum amorphous silicate 
feature. The black-body is aimed at contributing at wavelengths 
larger than 2Qfim and T^ont therefore constrained to be less than 
150 K. A typical value for the temperature T^om is found to be 
about 11 OK. 

In our fitting approach, only two free parameters charac- 
terize the continuum: the black-body temperature T^oni and a 
normalization offset (9v, at A2 = c/v2 ~ 13.5 yum. The imple- 
mentation of a variable offset comes from the fact that the 10 
and 18//m amorphous features partly overlap at A ~13-15yum 
(Fig-ID. and any realistic continuum should therefore pass be- 
low the observed flux at these wavelengths. For a given set of 
(^, Tcont and offset value (9v, , the synthetic continuum is obtained 
by solving the following two-equation system with respect to the 
normalization coefficients Kp\ and K\,\,: 



Kpi X A" + Khb X -Svi(rcont) - ^vi.obs 

Kpi X /Ij + Kbb X -Sv2(rcont) - -fvj.obs X (1 



Oy,) 



(2) 



where Fyobs is the observed spectrum, in units of Jansky. The 
constraints on the normalization coefficients {Kp\ > and K\,h > 
0) imply that for some objects, the best continua require K^b - 0, 
corresponding to continua described by a pure power-law. The 
normalization wavelengths are chosen to bracket the amorphous 
10//m silicate feature, with /1 1 = c/vi ~ 7.5yum, and/l2 - c/v2 ~ 
13.5/jm. 

The fit to the continuum-subtracted IRS spectrum is per- 
formed in two steps. First, a fit to the lOjum feature is obtained 
between Ai and A2 (warm component), then a second fit to the 
residual spectrum is obtained for wavelengths between Ai and 
/I3 ~35yum (cold component). The continuum-substracted IRS 
spectrum (/^v,obs - ^v.cont) around lOjum is reproduced within the 
range [Ai, At], by summing up the thirteen mass absorption co- 
efficients k! (A^species - 5 dust spccics, / indcx, and Assizes = 3 or 
2 grain sizes, j index depending on the lattice structure), multi- 
plied by a blackbody By{T^) at a given warm temperature r„, 
and weighted with relative masses Mi^ ■: 



Nsizcs Nspccics 

J=l '=1 



(3) 



The residuals (Fy.obs - ^v.cont - ^v.warm) are then similarly fitted 
between Ai and A^ with a synthetic cold component spectrum 
that writes: 



Nsizcs ^species 



^ ' size s " ■H^'^* 

i^v,coid = By(T,) X J] Yj "i^ Kr 



(4) 
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Fig. 2. Detailed example of B2C modeling steps for the object 
Sz 96. From top to bottom: the original spectrum and the esti- 
mated continuum (top panel), the continuum substracted spec- 
trum and the fit to the warm component (in orange, 2"'' panel), 
the residuals from the fit to the warm component and the fit to 
this cold component (in blue, 3"* panel), the final fit to the entire 
spectrum (in green) and the various contributions (4* panel), the 
relative residuals between the original spectrum and the final fit 
(last panel). 



The final synthetic spectrum, obtained in two steps, then reads: 



^ v,syiit ~ ^ v.coiit ' ^ v,waim "•" ^ i 



v.cold- 



(5) 



It depends on 2 x 13 parameters for the relative mass abundances 
of the warm and cold components (the M;^ . and M^ .), two tem- 
peratures {Tv/arm and Tcoid), and two parameters for the contin- 
uum (rcont and Ov,), which leads to 30 free parameters in total. 
The successive steps of the B2C modeling approach are de- 
composed in Figure|2]for the illustrative case of Sz 96. The orig- 
inal spectrum (Fv.obs) and the estimated continuum (Fv.cont) are 
displayed in the top panel. Then, the fit to the lOjum feature 
is performed (Fvwaim, orange curve on second panel) and sub- 
tracted from the continuum-subtracted spectrum (black line on 
the third panel). It shows that, although the contribution of the 
first fit is not negligible up to ~ 25 jum, a second, colder compo- 
nent is required to account for the large wavelength spectrum. 
The fit to the residuals (Fy^coid, in blue on the third panel) is 
then computed, leading to an overall fit (/^v,synt, green line) to 
the entire IRS spectrum as drawn in the fourth panel together 
whith all the contributions (continuum, warm and cold compo- 
nents). The last panel of Fig. |2] displays the final, relative resid- 
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uals for the full fit of the IRS spectrum, calculated as follows: 

l^v,obs ~ ^y,synt//^ v,synt- 



2.3. Fitting process 

The parameter space has a high dimensionality in our problem 
and imposes a specific fitting approach to appreciate the reliabil- 
ity of the results. We develop in this paper a method based on a 
Bayesian analysis, combined with a Monte Carlo Markov chain 
(MCMC)-like approach to explore the parameter space. 

Our procedure is built in order to randomly explore the space 
of free parameters (2x 13 dust compositions, 2 dust temperatures 
and 2 parameters for the continuum). We start with a randomly 
chosen initial set of parameters, then one of these parameters 
is randomly modified, while all the others remain unchanged. 
When the chosen parameter is related to the abundance of a grain 
species (the M^^ ■ and M^ .), the maximal value for the increment 
is 1% of the previous total mass for the considered component 
(warm or cold). When the chosen parameter is a temperature or 
the offset Ov2, we allow increments of at most 4% of the previ- 
ous temperature or the previous offset, respectively. These maxi- 
mum values were chosen to obtain small enough increments and 
therefore explore continuously the parameter space. 

For each set of parameters, a synthetic spectrum Fy^sym is cal- 
culated and the goodness of the fit to the observed spectrum is 
evalutated with a reduced x^- The parameter space exploration 
is therefore very much alike a MCMC approach, with two speci- 
ficities: we only perform jumps for oneparameter at a time, and 
second, the jumps are always acceptecU After « iterations (typ- 
ically n = 800), we set all the parameters to those that gave the 
lowest ;^'~ among the « previous iterations. This loop is done m 
times (usually m is also set to 800) and all the compositions, 
temperatures, offset and their associated ;if^ values are stored. 

The n x m reduced ;if^ values are transformed into probabili- 
ties assuming a Gaussian likelihood function (oc exp(-Y^/2)) for 
Bayesian analysis. Marginalized probability distribution for each 
free parameter are then obtained by projection of these probabil- 
ities onto each dimension of the parameter space. The best fit to 
the observed spectrum among all the simulations (i.e., the one 
with the lowest x^ value), yields relative masses for all the dust 
species and grain sizes {M^^^ and M^^^ parameters in Eqs.[3]and 
lU for both the warm and cold components, as well as best tem- 
peratures Tw and r^, which typically range between 200 K and 
300 K for the warm component, and between 70 A" and 150 A' 
for the cold one. The 1-cr uncertainties on the parameters are 
derived from the probability distributions for each parameter. 
We over-sample each probability distribution to compute half- 
width at half maximum for both sides of the distribution, and 
derive minimum and maximum 1-cr uncertainties. Figure[3]dis- 
plays ten examples of well -peaked probability distributions, with 
their respective over-sampled distributions overplotted, for some 
parameters of the fit to the Sz 96 spectrum, indicating the pa- 
rameters are constrained by the B2C fitting approach. Figure|4] 
displays four example spectra with their best fits (and residuals) 
over the entire IRS spectral range, showing that good fits can be 
obtained for spectra with different shapes. Spectral regions with 
high residuals mostly correspond to regions with low signal-to- 
noise ratio. 



' the implementation of a Metropolis-Hastings rule to decide on 
whether a jump should be accepted or rejected with a certain probabil- 
ity did not improve the fitting process, while increasing the calculation 
time 



The robustness of the B2C procedure has been intensively 
tested and this work is reported in AppendixlAl Using theoretical 
spectra we search for any possible deviations to the input dust 
mineralogy, that could either affect the inferred grain sizes or 
crystallinity fractions. The main result is that even if there is a 
slight deviation for a few individual cases, in a statistical point 
of view for a large sample, the B2C procedure is robust when 
determining the dust mineralogy. 

3. Spectral decomposition of 58 TTauri IRS spectra 
with the B2C model 

We run the B2C compositional fitting procedure on 58 differ- 
ent stars (most of them being TTauri stars, except BD+31 634 
which is an Herbig Ae star), for which we obtained Spitzer/IRS 
spectra as part of the c2d Legac y program. The spectra are pre- 
sented in lOlofsson et alj ( l2009l) and we refer to that paper for 
details about data reduction. T he selection of the 58 objects out 
of 96 in lOlofsson et al.l (|2009|) is based on several criteria. First, 
some objects do not have Short-Low data therefore the amor- 
phous 10/im feature is not complete. Second, as the goal is to 
determine the dust mineralogy, we do not run the procedure for 
objects that do not show clear silicate emission features, or with 
pecuUar spectra (e .g. cold disks like LkHa330 or CoKuTau/4, 
iBrown et al.l 120071) . Finally, objects for which continuum esti- 
mation in the 5-7.5 yum spectral region was not possible using a 
power-law were eliminated. 

The spectral range used for the fits is always limited to a 
maximum wavelength of 35 fim. The first reason of this choice 
is that the products of the c2d extraction pipeline are, for our 
sample, in average limited to 36.6jum. In addition to this, the 
degrading quality of the end of the spectra, likely caused by the 
quality of the relative spectral response function used at these 
wavelengths, may lead to an over-prediction of the crystalline 
content, for the cold component. In a few cases, the spectra are 
rising for wavelengths larger than 35 //m, and the fitting proce- 
dure tries to match this rise with crystalline features (which are 
the only features strong enough in this spectral range). Finally, 
longer wavelengths may probe an even colder dust content and 
this would require the implementation of a third dust component 
in order to reproduce the entire spectral range. For all these rea- 
sons, we choose to limit the modeling to wavelengths smaller 
than 35 yum. 

The source list and relative abundances for every object can 
be found in Table lA. l1 Because of some degeneracy between 
amorphous olivine and pyroxene opacities, we sum their respec- 
tive abundances to a single amorphous component. The final fits 
to the 58 IRS spectra are displayed in Figs lA.614A.8l 

3.1. Grain size properties 

3.1.1. Mean mass-averaged grain sizes 

With the outputs of the "B2C" procedure for the 58 objects pro- 
cessed we have statistical trends on typical grain sizes necessary 
to reproduce the spectra. The mean mass-averaged grain sizes 
for the warm and cold components, (awarm) and (Ocoid), respec- 
tively, are calculated as follows: 

-1 



^N^i 



N^warm/cold/ 



N™, 



/'N.izes Nsr 



Zfl V M-' X V V M-' 
J / 1 w/c,/ / 1 / I w/c,/ 



(6) 



where a\ - 0.1 //m (small grains), 02 = 1.5 //m (intermediate- 
sized grains) and 03 = 6;um (large grains), with the masses 
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Fig. 5. Mean grain sizes for the cold component as a function of 
the mean grain size for the warm component 



as defined in Eqs.[3] and HI We further define (a™° , ,,) and 

^ ' — ' ' — ' ^ warm/cold' 

^"^warm/coid^ ^^^ mass-avcragcd grains sizes for amorphous and 
crystalline grains, respectively, in the warm or cold component. 
We obtain a mean mass-averaged size of (flwaim) - 2.28 //m 
for the warm component, and a comparable value of {flcoid) - 
2.02jum for the cold component. Figure|5] shows that mass- 
averaged grain sizes for both components are uncorrected with 



each other In order to quantify the strength of the correlation, 
we compute the Kendall t correlation coefficient and its asso- 
ciated probability P. The t value denotes if there is any corre- 
lation or anti-correlation (t = 1 or -1, respectively), and the P 
probability corresponds to the significance probability (from 
to 1, from the most to the less significant). Regarding the latter 
two quantities we obtain a Kendall r value of -0.013 and a sig- 
nificance probability P = 0.888. This overall suggests that the 
warm and cold disk regions considered in this study are inde- 
pendent, as both components show uncorrelated grain sizes in 
the inner and outer regions. This suggests that differential grain 
growth is not the sole process explaining the observed variations 
fr om object to ob j ect. T his result is in line with the conclusions 
of lOlofsson et alj ( l2009b and supportin g the B2C model assu mp- 
tions. In their study of 65 TTauri stars. ISargentet al] ( l2009h find 
that grains are larger in the inner regions compared to outer re- 
gions, and argue this difference can be explained by faster grain 
coagulation in the inner regions, where dynamical timescales are 
shorter In our study, we find no strong evidence of such a differ- 
ence (2.28 versus 2.02//m). However Fig. |5] shows a larger dis- 
persion in grain sizes for the warm component compared to the 
cold component. This could be a consequence of shorter dynam- 
ical timescales in the inner regions where grains are not frozen 
and will be strongly submitted to both coagulation and fragmen- 
tation processes, compared to the outer regions. 

Additionally, we investigate the different mass-averaged 
grain sizes, for both the crystalline and amorphous grains. For 
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3.1 .2. Mass-averaged grain size versus disk flaring 
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Fig. 6. Flaring indices Ft,qIFit, as a function of the warm mass- 
averaged mean grain size. 



As in lOlofsson et alj ( |2009|) . we find a trend between (a„arm) 
and disk flaring as measured by the flux ratio F^o/Fiy, (fluxes 
in Jy integrated between 13 + 0.5 ;um for F^, and 30 + l/^m 
for Fio). As illustrated in Figure|6] we find that large grains 
in the warm component are mostly present in flat disks, while 
smaller grains can be found in both flared or flat disks. This anti- 
correlation for the warm component has a Kendall correlation 
coeflicient t of -0. 3 1 with a significance probability of 5 .Ox 1 0""*. 
This anti-correlation is still present, when considering only the 
warm amorphous grain sizes (r = -0.32 with P = 3.5 xlO~^). A 
similar anti-correlation has also been found by iBouwman et al.l 
(|2008|) . and by Wat son et al., (.2009a.) in their Taurus-Auriga as- 
sociation sample. 



the warm component, crystalline grains have a mean mass- 
averaged grain size of (flwann) - 1.14yum while amorphous 
grains have (a™°„) - 2.50//m. For the cold component, we find 
KZd> = 0.79 jum and (fl™°) = 2.40 ^/m. As we did not include 
large crystalline grains in the fitting process, it is not surprising 
to obtain smaller mass-averaged sizes for the crystalline grains 
compared to amorph ous grains. How e ver thi s trend is supported 
by the results from iBouwman et alj ( l2008h . for seven T Tauri 
stars. 



Considering {a^^^^} as a function of the flaring degree 
(FjqJFh), we search for a similar result as in lOlofsson et al.l 
(20091) where we found that small crystalline grains are pref- 
erentially seen in flattened disks, while large crystalline grains 
can be found in a variety of flat or flared disks. We did not find 
a similar trend from the outputs of the modeling, the correlation 
coefficient being t = -0.040, with a significance probability of 
0.66. According to this model, the size of the cold crystalline 
material seems to be strongly unrelated with the flaring degree 
of disks. 
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Fig. 7. Grain sizes distribution, for warm (red circles and dotted 
line) and cold (blue squares and dashed line) component. The 
MRN reference distribution (-3.5 index) is the black plain line. 



Fig. 8. Crystalline fraction of the cold component (Ccoid) as a 
function of the warm crystalline fraction (Cwarm)- 



3.1.3. Flattened grain size distributions 

Since our spectral decomposition includes three grain sizes, for 
the amorphous species, we can evaluate the shape of the size 
distribution in the atmospheres of disks. Assuming a differential 
size distribution d«(a) = Kaa''da where Ka is a normalization 
constant, we approximate the number Na of grains of size a with: 



Na - I d«(fl) 

Ja-Sa/2 

^ KaO^'da 



Kg 

l+P 



Sa\^^P 



da\^^P 



Ht) -("-?) 



(7) 
(8) 



after a Taylor expansion to the first order of Eq.Qabout a assum- 
ing 6a <s: a. This shows that NgJNat - {oilak)'' for two grain 
sizes fl, and a^, and therefore that relative numbers of grains can 
directly be used to evaluate the slope of the size distribution p. 
For each grain size, we compute the mean mass obtained from 
the B2C simulations, and then divide it by the corresponding 
volume (oc fl^). We then normalize these relative grain numbers 
so that the total number of 0. 1 //m-sized grain equals 1 . 

Fig. |7] shows the mean differential grain size distribution in 
normalized number of grains obtained following this procedure. 
On Fig.|2l the warm component is represented in filled circles, 
the cold component in open circles, and a reference MRN differ- 
ential size distribution {p - -3.5) with a dashed line. Assuming 
power-law size distributions, we find p indexes of pwarm =-2.89 
and /^coid =-3.13 for the warm and cold components, respec- 
tively, indicating much flatter size distributions compared to the 
MRN size distribution. Two additional runs of the B2C proce- 
dure for the 58 objects allowed to confirm this trend, and to es- 
timate uncertainties on the slopes. For the warm component, we 
typically find /9 = -2.90 + 0.1, and p = -3.15 +0.15 for the cold 
component. 

Because the p indexes are larger or close to -3, an immediate 
consequence of this result is that the emission in disk upper lay- 
ers is statistically dominated by the /jm-sized grains in our stellar 
sample, especially for the warm component, and is largely inde- 
pendent of the minimum grain size of the size distribution as 
long as its va lue is small enough ( submicronic, see discussion 
in Sec. 5.3 of lOlofsson et al.ll2009l) . This suggests that the flat, 
boxy lOyum feature profile for most TTauri stars discussed in 
terms of a depletion of small grains in Olofsson et al. (2009), is 
more precisely revealing a significant flattening of the size distri- 
bution, i.e. a relative lack of submicron-sized grains with respect 



to micron-sized grains, but not a complete depletion. This size 
distribution flattening is further discussed in Sec. ID 



3.2. Silicate crystals 

3.2.1. The crystallinity paradox reexamined 

Figure[8] shows the correlation between the warm and cold crys- 
talline fractions (which we denote as Cwaim and Ccoid, respec- 
tively, in the following), with Kendall's r =0.25 and a signifi- 
cance probability P = 5.8 xlO"-'. While slightly dispersed, there 
is a tendency for a simultaneous increase of the crystallinity in 
both the warm and cold components. The crystalline distribu- 
tions for the warm and cold components are displayed on the 
left panel of Fig.|9] For the warm component, the mean crys- 
talline fraction is (Cwarm) - 16%, while this fraction shifts up 
to (Ccoid) - 19% for the cold component. Overall, these results 
show no significant difference between the crystalline fractions 
in both components. 

These modeling results give us new insights concerning the 
crysta llinity paradox identified observationally in lOlofsson et al.l 
(|2009|) . The crystallinity paradox expresses the fact that crys- 
talline features at long wavelengths are ~3.5 times more fre- 
quently detected than those at shorter wavelengths. Using sim- 
ple models of dust opacities, we concluded that a contrast ef- 
fect (i.e. the strong lOyum amorphous feature masking smaller 
crystalline features) could not be accounted for the low detec- 
tion frequency of crystalline features in the lOfim range, a result 
valid for crystallinity fraction larger than 15% for 1.5yum grains 
(lOlofsson et al.ll2009l) . For lower crystallinity va lues, we indeed 
showed that the synthetic spectra computed in lOlofsson et al.l 
(2009J) were not representative of the observations. This contrast 
effect could therefore not be constrained for low-crystallinity 
fractions, where most of our objects lie for the warm component, 
according to Fig.|9] With the outputs of the modeling procedure, 
we can investigate this issue with a new look. Right panel of 
Fig. |9] shows the cumulative fractions as a function of the crys- 
tallinity fractions. Even if there are a few differences between the 
warm (red dotted) and cold (blue dashed) components, the two 
cumulative fractions display a very similar behavior that cannot 
solely explain the crystallinity paradox that we derived from the 
observations. This therefore means that, with this modeling pro- 
cedure, a contrast issue around lOjUm is requi red to match the 
observations described in lOlofsson et al.l(l2009l) . In other words. 
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Fig. 9. Left panel: crystalline distributions for warm component (Cwarm, red circles and dotted histogram) and cold component (Ccoid, 
blue squares and dashed hsitogram). Right panel: cumulative fractions of the crystallinity fractions, for the warm (red dotted line) 
and cold (blue dashed line) components. 



for some objects, few crystalline grains are required to repro- 
duce the 10 /im feature, but these grains do not produce strong 
features (on top of the amorphous feature) that can easily be de- 
tected in the spectra. Still, the left panel of Fig.|9]shows that the 
cold crystalline distribution is wider than the warm distribution. 
For the warm component, 3.5% of the objects have a crystallinity 
fractions larger than 40%, while this fraction shifts up to 13.8% 
for the cold component. This overall means that in a few cases, 
we see more crystalline grains in the cold compon ent compared 
to the warm component, as result also found by IS argent et al.l 
(l2009l) 



3.2.2. Crystallinity versus disk and stellar properties 

We also search for correlations between disk flaring proxies and 
crystallinity. We find no striking correlations regarding the disk 
flaring indexes Ft,o/F[t, and warm crystalline fraction Cwaim (t - 
-0.11 with a significance probability P = 0.23). Flared or flat 
disks present a wide range of crystalline fractions, with a strong 
dispersion. We find similar results, with an important dispersion, 
for the flaring index versus cold crystalline fraction Ccoid (t - 
-0.16 with a significance probability P = 0.08), meaning that 
crystallinity fraction of the cold component does not strongly 
depend on the disk flaring. 

Figure[TO]displays the dependence of crystalline fractions for 
both the warm (Cwaim, red open circles) and cold components 
(Ccoid, blue open squares) as a function of the spectral type, for 
stars between KO and M5. We find no correlation for the warm 
component with the spectral type (r = 0.13 with P - 0.20), sug- 
gesting that the degree of crystallization does not depend upon 
the spectral type (in the explored range), and that crystalliza- 
tion processes a re very general for T Ts. This result is in good 
agreement with IWats on et al.l (l2009ah . who find no correlation 
between crysta flinity and ste l lar lum inosity or stellar mass. On 
the other hand, IWatson et al.l (l2009al) only studied the cold disk 
regions crystallinity via the presence of the 33.6jum forsterite 
feature. Using the outputs of our B2C procedure, we are able 
to better quantify the crystalline fraction of the cold component, 
and we do not find any correlation of Ccoid with spectral type, 
the dispersion being too large (t - 0.12 with P - 0.24), thereby 
confirming the result bv lWatson et alJ (12009 ah . 
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Fig. 10. Crystalline fraction for the warm (Cwarm, open red cir- 
cles) and the cold (Ccoid, open blue squares) components, as a 
function of the spectral types. 



3.2.3. Dust evolution: coagulation versus crystallization 

Dust is evolving inside circumstellar disks, either in size or in lat- 
tice structure, in particular through coagulation and crystalliza- 
tion. To investigate possible links between these two phenom- 
ena, we search for correlations between grain sizes, for amor- 
phous and crystalline dust contents ((awarm/coid)), as a function 
of crystalline fractions, in both the warm and cold components 

v^warm/cold/- 

We find no evidence that grains are growing and crystallizing 
at the same time in the warm component. Considering (flwaim) 
and Cwarm values, we obtain the following coefficients: r = - 
0.038, with a significance probability equal to 0.67. This result 
is displayed in Fig.[TT] We obtain an even less favorable corre- 
lation coefficients for the cold component (t - 0.013 with P = 
0.89). This overall indicates that the processes that govern the 
mean size of the grains and their crystallinity are independent 
phenomena in disks. 

We also search for correlations between the amorphous and 
crystalline dust contents for the two components, to see how they 
evolve with respect to each others. Based on our B2C spectral 
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Table 1. Cloud-to-cloud variations for several parameters from 
the spectral decomposition. 



Fig. 11. Warm crystalline fraction Cwarm versus warm mass- 
averaged grain size (flwai-m)- 
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Fig. 12. Mean composition for warm and cold components ac- 
cording to the results of the 58 objects processed. 



decomposition, we find no significant correlations between the 
several amorphous or crystalline dust populations. 

3.2.4. Enstatite, forsterite and silica 



iBouwman et al.l ( |2008[) find that enstatite dominates in the in- 
nermost warm regions over forsterite, and on the other side that 
forsterite is the dominant silicate crystal in the cooler outer re- 
gions. Fig. [12] shows the mean composition we obtain for the 
58 objects, for both components (warm and cold). As amor- 
phous pyroxene and amorphous olivine spectral signatures are 
very similar, we merge together their relative contributions. 
Regarding the crystalline content, even if there is no clear pre- 
dominance of forsterite or enstatite masses for the warm com- 
ponent (8.7% vs. 7.4% respectively), for the cold component 
forsterite appears more abundant than enstatite (11. 4% vs. 7.4% 
in mass). In order to have more details on the presence or ab- 
sence of these two crystals, Figure[T3] shows the distribution for 
the enstatite and forsterite abundances (left and right panel, re- 
spectively), for the warm and cold components (filled circles and 
open circles, respectively). No predominance of either enstatite 
or forsterite is found in the warm component, while forsterite 
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seems to be more frequent than e nstatite in the cold c omponent. 
Similar distributions are found bv lSargent et alj (|2009|) regarding 
the warm and cold forsterite and enstatite. 

We search for correlations between the masses of the warm 
or cold enstatite populations, and warm or cold forsterite popu- 
lations. We find a correlation between warm enstatite and warm 
forsterite grains, with a Kendall's t value of 0.33 and with a sig- 
nificance probability P equal to 2.7x10"'*. This trend tends to 
confirm the results displayed in Fig. [121 the crystallisation pro- 
cesses in the warm component do not favor the formation of one 
or the other crystals. We also find a correlation between the rel- 
ative masses of warm enstatite grains and cold forsterite grains, 
with T = 0.31 and P equal to 7.0x10"'*. This tends to indicate 
that the enrichment in crystalline grains in the disk is global 
within ~10AU, and not local. 

Knowing that forsterite and silica can interact with each 
other to form enstatite (Mg2Si04 + SiOa -^ 2MgSi03), we 
search for correlations between the relatives masses of forsterite, 
enstatite and silica in both components. We may expect to find 
an anti-correlation between the presence of enstatite and silica, 
as silica disappears from the disk medium in the above reaction. 
We did not find any striking trend between all the masses, even 
when considering only small and intermediate-sized grains (0. 1 
and 1.5yL(m). This suggests that the reaction between forsterite 
and silica is not the main path for enstatite production. 



3.3. Sample homogeneity 

The studied stellar sample is drawn from several star forming 
regions. It may therefore be interesting to see if any cloud-to- 
cloud variations can be observed. The 58 Class 11 objects are 
distributed among 6 different clouds: Chamaleon, Ophiuchus, 
Lupus, Perseus, Serpens and Taurus (except for IRAS 08267- 
3336 which is an isolated star). Table[T]shows several results as 
a function of the corresponding clouds: the number of objects, 
the grain-size distribution slopes for the warm and cold compo- 
nents, and both the warm and cold crystallinity fractions. The 
results for the Perseus, Serpens and Taurus clouds are not sta- 
tistically significant, given the low-statistics number of stars for 
these regions. On the other side, the three other clouds show very 
similar behavior for all the considered results. The grain size dis- 
tributions and crystallinity fractions are very close to each other 
(with a number of objects between 13 and 16 per cloud). 

Recently, lOliveira et al.l (1201 Ol) studied the amorphous sili- 
cate features in a large sample of 147 sources in the Serpens, 
and compared their results with t he Taurus young star s. That 
study, along with the results from lOlofsson et akl (l2009l) . show 
that the 10 /im feature appears to have a similar distribution of 
shape and strength, regardless of the star forming region or stel- 
lar ages. In terms of silicate features, the stellar sample analyzed 
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Fig. 13. Distribution for the warm (red circles and dotted histogram), cold (blue squares and dashed histogram) for enstatite {left 
panel) and forsterite {right panel) crystalline abundances. 



in this study is therefore representative of a typical population of 
TTauri stars. 



4. Discussion 

4.1. Implications for the dust dynamics in the atmospheres of 
young disks 



In lOlofssonet al.l (l2009h . we discussed the impact of a size dis- 
tribution on the shape of the 10 //m feature, and could not disen- 
tangle between two possibilities to explain the flat, boxy 10 /im 
feature profile that is characteristic of micron-sized grains. One 
possibility was that the grain size distribution is close to a MRN- 
like distribution {p = -3.5), extending down to sizes of at least 
several micrometers, with the consequence that the Spitzer/IRS 
observations would then probe a truncation of the size distri- 
bution at the minimum grain size (of the order of a microme- 
ter). A second possibility implied a differential size distribution 
dn{a) CSC a''da that departs from a MRN-like distribution, being 
much flatter and with an index p close to, or larger than -3 to get 
the dust absorption/emission cross sections dominated by grains 
with size parameters Ina/A close to unity. 

By model fitting 58 IRS spectra of TTauri stars, we find a 
general flattening of the grain size distribution (pwarm - -2.9 
Pcoid - -3.15) in the atmospheres of disks around T Tauri stars 
with respect to the grain size distribution in the ISM up to about 
10 AU (this effect being stronger for the inner regions of disks). 
Our results suggest that the frequently observed boxy shape of 
the lOfim feature is rather due to the slope of the size distri- 
bution, than to a minimum grain size in the disk atmospheres 
close to a micrometer. This finding very much relaxes the need 
for a sharp, and quite generic truncation of the size distributio n 
at around one micrometer as proposed in lOlofsson et al.l (l2009h . 
Nevertheless, radiation pressure and/or stellar winds which had 
been proposed to explain such a cut-off' are likely to operate any- 
way, and may thus contribute to the flattening of the size distribu- 
tion by removing a fraction of the submicron-sized grains. This 
would furthermore help transporting crystalline grains outwards. 

An alternative explanation for the size distribution flatten- 
ing involves an overabundance of micrometer-sized grains with 
respect to the submicron-sized grains. At first glance, this situa- 
tion may sound counterintuitive as micrometer-sized grains are 
more prone to settle toward the disk midplane which would, in 



turn, tend to steepen the size distribution in disk atmospheres. 
This ignores that both turbulent diff'usion and collisions be- 
tween pebbles may complicate this picture, and may even pos- 
sibly revert this trend by supplying the disk atmospheres with 
micrometer-sized (and larger) grains. Therefore, the IRS obser- 
vations may in fact indicate that coagulation and vertical mixing 
are slightly more efficient processes to supply grains in disk at- 
mospheres than fragmentation, resulting in an overabundance of 
micrometer-sized grains compared to an ISM-like distribution. 
This situation may last over most of the TTauri phase. 

Interestingly, iLommen et al.l (l20I0h also found a flattening 
of the grain size distributions for T Tauri stars, probed by the 
millimeter slope of the SED. They found a tentative correlation 
between the shape of the 10 /im feature and the nun slope for a 
large sample of sources, with flatter 10 //m feature being asso- 
ciated with shallower mm slopes. Similarly, iRicci et al] (l20I0h 
also find grain size distributions flatter than in the ISM, with 
3 mm PdBI observations of T Tauri stars in the Taurus- Auriga 
star forming region. Therefore, the size distribution flattening 
in disk atmospheres may extend to deeper and more distant re- 
gions, and to larger grains sizes than those probed by mid-IR 
spectroscopy. 



4.2. Silicate crystallization and amorphization 

In this study, we have been able to quantify the f raction of crys- 
tals fo r both the warm and cold components. In lOIofsson et al.l 
(l2009h . we identify a crystallinity paradox expressing the fact 
that crystalline features are much more frequently detected at 
long wavelengths compared to short wavelengths. By building 
crystallinity distributions with the results from the spectral de- 
composition, we find that even if the cold crystalline distribution 
is slightly broader than the warm distribution, they both show a 
similar behavior (see right panel of Fig|9]l. Therefore, this may 
suggest that the observational crystallinity paradox could be 
the consequence of a combined effect: relatively low crystalline 
fractions (< 15-20%) for a great number of the TTauri stars, and 
a contrast issue at short wavelengths preventing the direct iden- 
tification of weak crystalline features on top of the lOyum amor- 
phous feature. Nevertheless, it is noteworthy that correcting for 
the contrast issue through B2C spectral decomposition does not 
revert the trend, but instead makes the distribution of crystals in 
the warm and cold disk regions quite similar 
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In order to explain the non negligible crystallinity fractions 
obtained for the cold component (up to 70% in a few cases, e.g. 
Fig.lSJ, there is a strong need for radial outward transport, as- 
suming that the only source of production of crystals is thermal 
annealing or gas-phase condensation/annealing, taking place in 
the inner regions of disks. Several models have shown that such 
mechanism is likely to happen in circumstellar disks around ac- 
ti ve stars (see Keller & Ga il 2004 or Ciesla 2009). According 
to ICieslal ( |2009|) . for an accretion rate larger than 10^'' M^/yr, 
grains can be transported to distances out to 20 AU in about 
~ 10^ years. A res ult in line with the 2D time -dependent disk 
model described in lVisser & Dullemondl(l2010l) . However given 
the relatively young ages of our stellar sample (~ few Myrs), and 
the rather high accretion rate required (10^^ M^/yr), we cannot 
assess that radial transport is the only mechanism responsible 
for the strong similarity of both warm and cold crystallinity dis- 
tributions. This instead suggests that other crystallization pro- 
cesses are necessary, that must take plac e in the outer region s 
of disks. A theoretical model, detailed in iTanaka et al.l (1201 Oh . 
shows that depending on the gas density, crystallization can take 
place with typical temperature of a few hundred Kelvin. Such 
process can be triggered i n the outer regions by nebular shocks 
dPesch & Connollvll2002l) which lead to graphitization of the 
carbon mantle providing the latent heat for silicate crystalliza- 
tion. 

However, crystalUzation by thermal annealing is likely to 
happen on very short timescales (few hours at T ~ lOOOK, 
according to laboratory measurements, see iHallenbeck et al.l 
Il998h . Therefore, one can expect to observe higher crystallinity 
fractions in the inner warm regions compared to the outer cold 
regions. One possible explanation to this issue may rely on the 
importance of accretion activity and on amorphization processes 
of crystals that are preferentially tak ing place in the in ner regions 
of d isks, as tentatively pro posed bv iKessler-Silacci e t al. (2006) 
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Fig. 14. Crystalline fraction for the warm component (Cwarm) as 
a function of the X-ray luminosities. 



l2009bh with the X-ray emission from one or two-temperature 
plasma combined with photoelectric absorption, for a subset of 
17 objects among our stellar sample analyzed with the B2C com- 
positional fitting model. The low overlap between X-ray obser- 
vations with the c2d sample can be explained by the different 
observation strategies; X-ray campaigns are mostly targetting 
dense cores while our observations are located at the periph- 
ery of the clouds. Fig.[T4]displays the crystalline fraction Cwarm 
for the warm component inferred with the B2C procedure, as a 
function of the X-ray luminosity. The associated Kendall t value 
is 0.16 with a significance probability P - 0.36, suggesting a 
dispersed weak trend, if any. To confirm any weak trend, we 

^ would need to have more measurements below about 2 x 10^^ 

and Olofsson et al.l (1200 ^. and recently discussed in more de- erg.s"'. Above this luminosity value no trend is visible. Still, we 

observe that objects with high X-ray luminosities seem to have 
higher crystalline fracti ons, which seems at fi rst in contradiction 
with the conclusion bv lGlauser et al.l (l2009l) . Given the limited 
number of objects for which we have X-ray data in our sample, 
however, we were n ot able to perform any age selection as in 
fClauseret all ( 120091) . 



tail in Glauser et al.l (l2009i) . This possibility is further supported 
by the results of lAbraham et al.l (|2009|) who studied the TTauri 
star EXLup with IRS, at two different epochs separated by a 
three-year time span , first i n a quiescent phase (see spectrum in 
iKessler-Silacci et al] l2006h . and then two months after an op- 
tical outburst driven by an increase of the mass accretion rate. 
They find an increase of crystallinity caused by this outburst in 
the inner disk regions, supposedly due to thermal annealing of 
amorphous silicates. The authors did not find any colder crys- 
tals in the outer component after the outburst, and therefore ex- 
cluded shock heating at larger distances. This study shows that, 
obviously, the accretion activity may have a stro ng impact on 
the du st chemistry in circumstellar disks. Recentlv. lGlauser et al.l 
(|2009|) identified an anti-correlation between X-ray luminosity 
of the stars and crystalline fraction for the warm disk component 
(Cwarm in Sec l3.2b . The authors performed a decomposition fit 
to the lOyL/m feature for a stellar sample of 42 TTauri stars, and 
find an anti-correlation between Cwaim and X-ray luminosities 
for a subsample of 20 objects with ages between 1 and 4.5 Myr 
This anti-correlation suggests that the stellar activity and ener- 
getic ions from stellar winds, traced by X-ray flux, can amor- 
phize crystalline grains in disk atmospheres due to high-energy 

particles. 

In a similar approach as iGlauser et al.l (l2009l) . we attempt 
to investigate the impact of high-energy particles, assumed to 
be traced by X-ray luminosity, on the disk mineralogy. Table|2| 
displays the X-ray luminosities (in units of erg. s"') available in 
the literature or obtained by fitting the spectra fro m the Second 
XMM-Newton serendipitous source catalogue (IWatson et al.l 



4.3. Strength of the Wfim feature, a proxy for gram size 7 
The strength ^p^^™ of the 10 /^m feature has been com- 



^ reaK i ■ i ■ 

monly used in pr evious studies ('e.g. IvanBoekel et al.l I2005L 
IKessler-Silacci et al. 2006, Olofsson et al.ll2009l) as a proxy for 
grain size, small S^^^ values tracing large (micrometer-sized) 

grains, while larger S^^"^ woul d be tracing s r naller ( submicron- 
sized) grains. In a recent paper. IWatson et aLJ ( l2009al) instead at- 
tribute the diversity of strengths of the 10 /im feature to a com- 
bination of sedimentation and an increase of the crystallinity de- 
gree (as in Sargent et al. 2009). With the outputs of our B2C 
spectral decompositional procedure, we revisit the interpretation 
of the strength of the 10 /im feature as a proxy for grain size. 

The left panel of Fig.[T5]displays the correlation between the 
mass-averaged grain size for the warm grains ((flwaim)) and the 
strength of the lO/^m feature (•S'pj,^™). We find a Kendall t cor- 
relation coefficient of about -0.30 with a significance P prob- 
ability smaller than 1.1 xlO""', showing that grain size is well 
trac ed by the strength of the feature. It is noteworthy to point out 
that ISargent et alj ( l2009l) find a different result regarding these 
two quantities, in their spectral decomposition analysis of 65 
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Table 2. X-ray luminosities for a subset of our stellar sample. 
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Star Name 


Satellite 


Lx [erg/s] 


Ref 


HM27 


ROSAT 


< 6.3x10-'' 


Feiselsonetal. (1993J 


LkHA271 


CXO 

XMM 


1.3x10^'* 
1.6x1029 


Getman et al. (2002) 


Hn9 


Robrade & Schmitt (2007) 


TWCha 


ROSAT 


4.0x10^' 


Feieelsonetal.n993) 


IMLup 


XMM 


4.2x1029 


this study 


CK4 


CXO 


4.4x1029 


Giardino et al. (2007) 


SSTc2dJ162715.1-245139 


XMM 


4.6x1029 


Ozawa et al. (2005) 


SXCha 


XMM 


4.7x1029 


this study 


XXCha 


XMM 


l.OxlO'o 


Robrade & Schmitt (2007) 


VZCha 


XMM 


1.3x10^" 




V710Tau 


XMM 
XMM 


1.4x10^° 
1.5x1030 


GUdel et al. (2007) 


WXCha 


Stelzer et al. (2004) 


SR9 


ROSAT 


2.5x10^° 


Casanova etal. (1995) 


RULup 


XMM 


3.1x10^" 


this study 


C7-11 


XMM 


5.9x1030 


Stelzer et al. (2004) 


Sz96 


XMM 


6.2x1030 


this study 


VWCha 


XMM 


9.3x1030 


Stelzer et al. (2004) 
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Fig. 15. Left panel: Mass-averaged grain sizes for the warm component, as a function of the strength of the lOyum feature (SpJ^'^) 
Right panel: power-law index p of the grain size distribution as a function of S „ ^'/'. 



T Tauri stars. As show n in the review by IWatsonI ( l2009l) , based 
on lSargent et al.l(l2009f) results, they find a correlation coefficient 
of 0.08 with a significant probability of 0.55. This means, ac- 
cording to their model, that the strength of the 10 //m feature is 
not correlated with the presence of large grains. More precisely, 
they find that seven objects among their sample with the small- 
est lOyum feature strength, can be reproduced with low abun- 
dances of large grains. However, our opposite result is supported 
by several points: we successfully checked for all dependencies 
that could lead to a possible overestimate of the grain sizes (e.g 
the continuum offset (9v,) and we show in App. lA.ll based on 
the analysis of synthetic spectra, that we may eventually under- 
estimate the grain sizes in the warm component rather than over- 
estimate it (underestimated by 7 + 1 1 % with respect to the input 
value). This deviation ((fl)/(flinp)) is furthermore independent of 

the strength of the lO/^m feature (S p^^"), with a r value of -0.06 
and a significance probability of 0.52. 

We also search for a correlation between the warm crys- 
talline fraction Cwarm and the strength of the 10 /im feature 
(S p '^f'), and find a highly dispersed relation between the two 



quantities. Fig. [T6l shows that the distribution is flat (r - -0.15 
with a significance probability of 0.10), which suggests that 
crystallization cannot be the sole explanation for the diversity 
of lOyum feature stren gths. This is in line with the results by 
iBouwman et al.l (l2008h who do not find any coiTelation between 
these two quantities in their seven T Tauri star sample. 

As mentioned in lOlofsson et al.l (|2009|) and in this paper, the 
grain size distribution in the upper layers of disks has a strong 
impact on the 10/im feature strength. The right panel of Fig. [15] 
shows the power-law indexes p of the size distributions (cal- 
culated as in Sec. l3.1.3l l. as a function of the SpJ^'^ values. We 
find an anti-correlation (t - -0.28 with P = 2.3 xlO"^), indi- 
cating the shape of the 10 //m feature is indeed related to the 
slope of the grain size distribution. In fact we see two different 
regimes appearing: the small S pj^^^ values mostly correspond to 
objects with slopes p larger than -3 (dashed horizontal line on 
the figure), while objects with \aigeS pj^^^ values essentially have 
slopes p smaller than -3. This points in the same direction as 
previously discussed, namely that the variation in S pj^^ may to 
a large extent trace the variation in the slope of the size distribu- 
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Fig. 16. Crystalline fraction for the warm component (Cwarm), as 
a function of the strength of the lO^um feature (S „ '"^'). 



tion. Not surprisingly, we find a very strong correlation between 
the mean mass-averaged grain size and the slope of the grain size 
distribution, with t = 0.79 and a significance probability below 
10-38. 



5. Summary and conclusion 

In this paper, we present a large statistical study of the dust min- 
eralogy in 58 proto-planetary (Class 11) disks around young solar 
analogs. We develop a robust routine that reproduces IRS spec- 
tra over their entire spectral range (5-35 //m), with a MCMC- 
like approach coupled with a Bayesian inference method. This 
procedure explores randomly the parameter space and repro- 
duces IRS spectra using two independent dust populations: a 
warm and a cold component, arising from inner and outer re- 
gions from disks, respectively. This 2-component approach is 
support ed by our previous stat istical analysis of the silicate fea- 
tures in lOlofsson et al.l (l2009h . This B2C procedure has been 
tested over synthetic spectra giving, statistically speaking, coher- 
ent and reproducible results. In its current form, the B2C model 
allows to derive relative abundances for 5 different dust species, 
3 different grain sizes for the amorphous grains, 2 sizes for the 
crystalline grains, for both inner regions (< 1 AU) and outer re- 
gions (~ 10 AU) of disks. Based on the modeling of the 58 IRS 
spectra of TTauris stars, we find the following results: 

1 . The grain sizes in the inner and outer regions of disks are 
found to be uncorrected, reinforcing the idea of two inde- 
pendent components probed by IRS. However, crystallinity 
fractions in warm and cold components are correlated to 
each other, indicative of a simultaneous enrichment of crys- 
tals within the first 10 AU of disks. This suggests that dynam- 
ical processes affecting the size distribution are essentially 
local rather than global, while crystallization may rather 
trace larger scale phenomena (e.g. radial diffusion). Finally, 
we find that grain size and crystallinity fraction are two in- 
dependent quantities with respect to each other in both the 
warm and cold components. 

2. We quantify the growth of grains compared to the ISM. We 
find a significant size distribution flattening in the upper at- 
mospheres of disks probed by IRS, compared to the MRN 
grain size distribution. While the MRN size distribution has 



a power law index p - -3.5, we find flatter distributions 
for both the inner regions {{p} - -2.90 ± 0.1), and the outer 
regions {{p} - -3.15+0.15). This explains why the IR emis- 
sion of T Tauri stars studied using Spitzer-IRS is dominated 
by yum-sized grains, despite the possible presence of signif- 
icant amounts of submicron-sized grains. The S J^^ value 
turns out to be a proxy of the slope of the grain size dis- 
tribution. This finding, combined with recent similar results 
from observations at millimeter wavelengths, suggests that 
the size distribution flattening is not confined to the disk up- 
per layers but may extend to deeper and more distance disk 
regions, and larger grains. 

3. We reexamined the crystallinity paradox identified in 
lOlofsson et alJ (l2009h . by building crystallinity distributions 
for the warm and cold components. The mean crystallinity 
fractions are 16 and 19% for the warm and cold components, 
respectively. Even if the crystalline distribution for the cold 
component is slightly wider compared to the warm compo- 
nent, both distributions show a very similar behavior, sug- 
gesting the silicate crystals are well mixed within the first 
10 AU of disks around T Tauri stars. According to the B2C 
model, the 3.5 times more frequently detected crystalline 
features at long (20-30yum) than at short (~10yum) wave- 
lengths arises from the combination of rather low crystalline 
fractions (<15-20%) for many T Tauri stars, and a contrast 
effect that makes the crystalline features more difficult to 
identify at around lOyum where the strong amorphous fea- 
ture can hide smaller crystalline features. 

4. We see a trend where flared disks (high F^o/Fit, values) tend 
to preferentially show small warm grains while flat disks 
(small F^o/Fiy,) show a large diversity of grain sizes. We do 
not find any correlations between the disk flaring and the 
crystallinity, meaning that the presence of such crystals is 
not related to the shape of the disks. 

5. We find no predominance of any crystal species in the warm 
component, while forsterite seems to be more frequent com- 
pared to enstatite in the cold component. Regarding the dif- 
ferent dust compositions, we find no evidence of correlations 
between them, except for the warm enstatite and forsterite 
grains, suggesting that the crystallisation processes in the in- 
ner regions of disks do not favor any of the two crystals. We 
find no link between silica and enstatite as one may expect if 
the main path for enstatite production would be the reaction 
between silica plus forsterite. 

6. We do not find any striking correlation between spectral type 
and warm or cold crystallinity fraction. We do not either find 
any correlatio n between X- r ay activ ity and crystallinity. As 
suggested by IWatson et al.l ( l2009al) . the time variability in 
young objects may be responsible for erasing most of the 
correlations regarding dust mineralogy. 
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Appendix A: Validation 

A.1. Procedure validation 

We evaluate the robustness of our procedure by fitting syn- 
thetic spectra with known abundances, temperatures and con- 
tinua. This allows to check whether these quantities could be 
recovered when the spectra are processed with the B2C proce- 
dure. To this end, we used the best fit, synthetic spectra for the 
58 objects analyzed in this paper (and presented in Sec.|3]l, as 
fake, but representative observations for serving as inputs to our 
B2C procedure. We reconstructed these synthetic spectra using 
the outputs of the fitting process on the original data. Namely, the 
same continuum emission, the same relative abundances for the 
dust species and their corresponding temperatures. The uncer- 
tainties chosen for the synthetic spectra are those of the original 
observed spectra. 

Fig. lA. II shows the ratio between the output and input crys- 
tallinity, as a function of the ratio between the output and in- 
put grain sizes, for both the warm and cold components. The 
plot shows ratios of about 1 on both axis, especially for the cold 
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Fig. A.l. Results for the fits to the 58 synthetic spectra to test our 
"B2C" procedure. The ,Y-axis shows the dispersion in grain size 
(ratio of inferred over input mean mass-averaged grain sizes) and 
the y-axis the ratio between inferred and input crystallinity. Red 
open circles correspond to the warm component, and blue open 
squares to the cold component. The filled circles are mean val- 
ues, with error bars corresponding to standard deviations. 



component (blue squares on the figure), suggesting some best 
fits obtained with the procedure may not be unique. But from a 
statistical point view, the B2C procedure produces reproducible 
results as the mean ratios between output and input quantities 
remain close to 1 (vertical and horizontal bars on the figure). 

More precisely, the warm component crystallinity is slightly 
overestimated by 22% with a dispersion (i.e. the standard devi- 
ation) about this value of 30 %. The inferred cold component 
crystallinity is satisfactorily reproduced at the 2% level with re- 
spect to the input value, and the dispersion, 35%, is rather similar 
to that for the crystallinity of the warm grains. The mean mass- 
averaged size of the warm grains is slightly underestimated by 
7+11% with respect to the input value, while the mean mass- 
averaged size of the cold grains is well reproduced but with a 
larger dispersion ( 1±32%). 

A closer look at the dispersion of the calculated crystalline 
fraction as a function of the input crystallinity is shown in 
Figure lA!2] for the 58 fits to synthetic spectra, for both the warm 
(red open circles) and cold (blue open squares) components. 
The over-estimation of the warm component crystalline frac- 
tions seen in Fig. lA. II is mostly visible for objects with low- 
crystalline fractions (below 20%). For the cold component, it 
actually tends to be rather underestimated for large crystalline 
fractions (above ~40%). The overall trend for both the warm and 
cold components is that the dispersion raises with the decreas- 
ing crystallinity. This finds an explanation in the fact that objects 
showing high crystalline fractions display strong, high-contrast 
crystalline features that are less ambiguously matched by theo- 
retical opacities than objects with lower crystalline fractions. 

A.2. Influence of the continuum on the cold component 

The continuum estimation is a critical step when modeling spec- 
tra, especially for the cold component, and may contribute to the 
larger uncertainties on the infeiTed cold component crystallini- 
ties and sizes discussed in previous subsection. The continuum 
may affect the model outputs in the following way: a high-flux 
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Fig. A.2. Ratios between output over input crystallinity fractions 
as a function of the input crystallinity fractions, for the fits to the 
58 synthetic spectra. Red open circles are results for the warm 
component and blue open squares are for the cold component. 
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continuum in the 20-30 //m spectral range leaves very little flux 
to be fitted under the spectrum, possibly leading to a composition 
with few large and/or featureless grains. A low-flux continuum 
could, on the other hand, lean toward large/featureless grains to 
fill the flux left. We therefore examine if (whether) such trends 
are present in the results of our B2C procedure (or not) for the 
cold component. 

We quantify the level of continuum by integrating the 
continuum-substracted spectra between 22jjm and 35 /im {x- 
axis on Fig. lA.3l l. Low-flux continua have high integrated fluxes 
and are located on the right side of the plot in Fig. lA.3l the high- 
flux continua being on the left side. The left panel of Fig. lA.3l 
shows the crystalline fraction for the cold component, as a func- 
tion of the integrated flux left once the continuum is substracted. 
The trend low continuum - low crystalline fraction is visible for 
the highest integrated flux values. But the large dispersion for 
low integrated flux values indicates that no strong bias is in- 
troduced as we obtain cases with very high-flux continua but 
still with very low crystalline fractions. It remains that for the 
low-flux continua cases, we cannot possibly obtain very high 
crystalline fractions as a lot of flux needs to be filled to match 
the spectra. Indeed, amorphous grains are the best choice to ful- 
fill this requirement, therefore diminishing the cold component 
crystalline fraction. 

The influence of the continuum on the estimated mass- 
averaged grain size is shown on the right panel of Fig. lA.3l A 
very weak and dispersed anti-correlation is found between the 
two parameters. We obtain a t value of -0. 12 with a significance 
probability P - 0.17, meaning that the adopted shape for the 
continuum is not strongly influencing the inferred mean grain 
size. The trend is mostly caused by a few low-flux continua ob- 
jects modeled with large grains. 

To conclude, the continuum estimation is a challenging prob- 
lem for the cold component and its adopted shape will always 
have an impact on the results for the crystallinity and grain size 
at the same time. Still, we have checked that, statistically speak- 
ing, we are not introducing a strong and systematic bias with our 
simple, two free parameter continuum. 



A.3. Importance of silica and necessity for large grains 

Usually, Mg-rich silicates ar e considered for the dust m ineralogy 
in proto-planet a ry disk s (e.g^ Henning & Meeusl2009l) . But both 
lOlofsson et"an (l2009h and ISargent et all (120091) attribute some 
features in IRS spectra of young stars to silica (composition 
Si02). To gauge the importance of silica in our B2C composi- 
tional approach, we run the B2C model with and without silica. 
Many fits were improved adding silica in the dust population (as 
shown for one example in Fig. IA.5l the 20-22 //m range being 
the spectral range where the improvement is the more notice- 
able). As can be seen on the right panel of Fig. lA.4l showing 
the reduced ;if J for simulations with and without silica, many fits 
are improved when silica is included, demonstrating the non- 
negligible importance of silica for our B2C model. 

We have also critically examined the need for large, amor- 
phous, 6.0//m-sized grains in the B2C model as they are almost 
featureless contrary to the 0. 1 and 1 .5 jum-sized grains. We there- 
fore run the B2C procedure on the 58 same objects with only 
0.1 and 1.5yum-sized grains. Left panel of Fig. lA.4l shows that 
the mean value of the reduced x^ for all the simulations is aug- 
mented by 2.50 when using only two grain sizes, showing that 
the majority of the fits are improved using three grain sizes in- 
stead of two. 

One can also worry about the influence of the offset (9v,2 val- 
ues on the inferred grain sizes. As a large offset value leads to 
larger flux below the continuum-substracted spectrum, it may 
indeed favor larger grains. We therefore examined if there was 
any correlation between the fraction of large 6.0/im-sized grains 
and the Ov,2 values, and we quantified these relations with cor- 
relation coefficients. Considering the warm large grains and the 
(9v_2 offsets, we find a t value of 0.04 with a significance proba- 
bility of P = 0.69. For the cold large grain fraction and the Oy^i 
values, we obtain t - -0.03 and P - 0.70. This means that there 
is no significant influence of the offsets (9v,2 on the inferred grain 
sizes. 
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Table A.l. Dust composition derived using the "B2C" procedure for 58 objects 
from our sample. For each star, first line corresponds to the warm component 
abundances and second line corresponds to the cold component abundances. 
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Table A.l. continued. 
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15 9+*^ 
59^4+10.2 



TT" 

5 7+--*' 
9 3+2:3 
QQ+1.6 

+7.8 
2.7 
+3.1 
3.3 
QQ+1.6 
"■"-0.0 

T T+3.8 

•'■•'-1.8 

11.8!|^ 

7 -j+bA 

4 5+2-^ 
7.6!" 

^-0.0 
0+'-' 

+0.7 
0.0 

'-'■'-'-0.0 

H.4 

-1.2 
^+1.3 

-0.9 
3+2.4 

-2.1 
;+4.7 

-2.3 
17 ^+4.8 

1+'*-^ 

-M 

-M 
^-0.0 

H.9 
^-1.8 

16 1 +■*■'' 

oo+'-3 

'^■'^-0.0 

8 9+5'' 

-2.0 

4 2+2-5 

'-'■•'-0.0 
2 1+2-'' 

1 7+5-' 
^ ■ ' -0.0 

4 1+3-' 
2.3+?-^ 



'0 

-0^0 
2.9 



0+"-^* 

4:7+15 
0+'-^ 

"■'-'-0.0 

0+2-'' 
'-'■'-'-0.0 

0+5-3 
"■'-'-0.0 

9 Q+1.4 

^•^-2.4 

0+2-3 
'-'■'-'-0.0 

11.9+^J 
0+3-^ 

0+3-2 
"■"-Q.Q 

0+°-' 

0+2" 
-Q.Q 

"■"-0.0 

Q 7 + 1.0 

"■'-0.0 
11 1+3-5 

"■"-0.0 

"■"-0.0 
0+'-3 

"•"-0.0 

l + '2-'' 
"■^-0.0 

1+2-'^ 
"■^-0.0 

17 6+35-8 

i ' •"-0.0 

0+2-^ 

"■"-Q.Q 

7+''-5 

"■'-0.0 
"■"-0.0 

7 3+3-3 

'■•^-4.8 

0+'^-'^ 

"■"-Q.Q 
"•"-Q.Q 

0+2-' 

"■"-0.0 

0+'-' 

"■"-0.0 

0+'-° 
"■"-0.0 

1 Q+4.9 
^■^-0.0 

"■"-0.0 

1 9+4.4 

-O.Q 

2 1+2-^ 

^•^-1.2 

7 0+3-** 

' ^"-4.4 

0.0+2-;^ 



133^ 

QQ+l.S 

2 1+'-° 
^■^-0.0 

2.8+2-0 



70 



15.6+2-5 
18.1!^.o 

T T + 1.6 

■'■■'-0.0 

15.0!*-^ 

'■-'-3.2 
ll-S+ll-i 



IRS60 

ISO-ChaII54 

LkHA271 

LkHA326 

LkHA327 

RNO90 

ROX42C 

ROX43A 

RULup 

RXJ 130 1.0-7654 

RXJ1615. 3-3255 

RYLup 

SR9 

SSTc2d J033036.0+303024 

SSTc2dJ033037.0+303128 

SSTc2dJ033052.5-H305418 

SSTc2dJ033241.7-H3 11046 

SSTc2dJ161159.8-382338 

SSTc2dJ162148.5-234027 



0.0 
20.0 



+-^9'* 



17.7 



0.0 

+3.9 

3.5 

+11.4 

'57 t: + lt.6 

-"■-'-n.o 

0+3-3 
"■"-0.0 

7.4+4-' 



2.9 

+3.3 
-0.0 
+4.8 



+l:f 



7.0 
+3.6 

5.3 
+5.9 
-2.7 
+1.6 

+i:l 



0.0: 

15. 1_ 
12.1+: 
22 7 

20.5+4-4 

14.9+7:1 

29.7 
19.1 
11.3 

39-4_5.i 

30.1^^ 

q 1 +9.4 

-4^5 

15.0+3-^ 

7+9.0 

57.4^4-^ 
28.6!!^ 
57.8^4-2 

17 2+3:^ 
14 0+3" 

26 7+"'-° 

1 0+2.1 
^■'^-o.o 

29 5+^-5 

9 T + 1.4 

^■■'-0.0 

57 3+'5-8 

-"■•'-14.0 
T 9+0.9 

■'■^-o.p 

59 8+'-3-8 
•'^■''-11.6 

4gQ+11.8 
^"■"-9.0 
9 /1+I8.6 
^■^-0.7 
1 9+4.8 



1.9 
+2.0 
0.0 
3.0^3-2 

0+"^:^ 

"■"-0.0 
9 C + I.4 

^■'^-o.o 

99 c + 12.9 
^"^■-'-11.2 

o.o!p 

1 s+0.6 

8.0!2-| 

A+4.6 



8.0-, 

4.8+2-^ 

17.4 

6.0 

1.3 

1.3 



0.0 

+4.7 

2.7 

+1.2 

+(y 
-0.0 
+1.4 

9^5!!! 

0+3-3 
"■"-0.0 

7+'-^ 

14+8:^ 
^■^-0.0 

2 3+'-3 

z.j_j4 
"■^-0.0 

4 4+1.3 
^■^-1.2 
14 5+''-4 

^■"-0.0 

Q 4+1.0 

"■^-0.0 

1 6+°-8 
^■"-0.0 

3+''-3 
"■■"•-o.o 

1 9+0.8 
^■^-0.0 

Q j+10.6 
"■^-0.0 

28.0+^4^ 

2 2+'-^ 

^■^-0.0 

9 7+3.4 
^^ -0.0 



0.2+ 
5.3+ 
2.V 



0.0 
0.0 

H.2 

-0.0 

j+3.5 

"-0.0 

4 0+10.2 

^■"-0.0 

0.4!'-' 

0+^:^ 
"■"-0.0 

"■•'-0.0 
9 q+i.O 
^■^-1.7 

o^9!i^ 
^■4!l^ 

0+'-^ 

"■"-0.0 
11 9+4-4 

0+°-'' 

"■"-0.0 
0+1-8 
"■"-0.0 

3+°-5 

"■•'-Q.Q 

"■"-0.0 

0+2-4 
"■"-0.0 

0.0+5-2 



0.0^ 
0.0^ 



0.0 
+0.8 

0.0 
+0.8 

0.0 
+2.4 



0.0 



0.0 

17 Q+10.3 
^'•^-7.2 

i.9!!:2 

19.5!^-5 

1+'-^ 

"■^-0.0 

1 7+2° 
■ -QS 

"■"-Q.Q 

5 7+4- '^ 
-'■'-o.o 

13 4+3-3 

0+2-2 
"■"-0.0 

0+'-2 

"■"-0.0 
4 9+3.8 
^■^-3.4 

7+°-5 

"■ ' -Q.O 

"■"-0.0 

"■"-0.0 
1 + 1.8 

"■^-n.n 

0+'^-'^ 
"■"-0.0 

o 4+18.4 
"■^-0.0 

0+'-5 

"■"-Q.Q 

0+°-" 

"■"-0.0 

0+°-4 

"■"-0.0 

11 7+3-0 

i i . / 4 5 

2 7+1.7 
^■'-1.1 
0+'-4 

"■"-Q.O 
"■"-0.0 
"■"-0.0 
"■"-0.0 

5+''-^ 
-Q.Q 

"■"-0.0 

0+'3-5 
"■"-0.0 

^■(>-lt 

'■"SI 

-0.0 



ri33''2 



5.4 
22.5 
63.9 



0.0+ 



11.7 
+7.6 

1.3 
+8.5 

-4.0 

16.2 
16.6 

55.0+f,^ 
33.7+3:4 
28.3+4:4 
6.7^4-^ 

40.0^+4-^ 
23.5^3:1 

140+2:1 
^^■"-1.9 

40.0!^-2 
23.6+^8.2 

33.0^4-^ 

31 7+2^-' 
•' ^ ■ ' -0.0 

28.5!^-2 
45.1!3:^ 
20.2+^-^ 
28.6!g 
49.5^5-2 
24.9+5:3 
14.8!j; 

41 2+''-^ 
^^■^-7.2 

38.2+4-8 

23.8!« 
26.7!2.2 
99 9+8.5 

^^■^-0.0 

21 5+4-'' 

3 0+''-° 
•'■"-1.7 

H4.2 
-20.3 
+3.9 
-4.8 



2.0 



14.5 
7.7 



0.0+ 



0.0 



7.6+ 
6.0! 
3.9+ 

4.5+ 



0.0+ 

0.0! 

3.3+ 



6.8+ 

l.V 



2.3 
12.9 



0+" 
"■"-0.0 

+5.7 
12.0 
+2.1 

9-9 

+4.4 



45.7 

0.0 

0.0_ 

2 9+^:'^ 
^■^-2.3 

o ^+3.8 

'^■•'-4.0 

0+4-2 
"■"-Q.Q 

9 6+5-3 

^■"-2.2 

3.4!3-'' 
2 4+^:^ 

^■^-0.0 

10.3+4-^ 

10.7!f-3 
4 4+1.4 

10 9+'-5 
^"■^-1.0 

1 5+'-5 
2.0 

16.1_g3 

19.8+4-^ 
4 8+13.« 

Oo+8:§ 
"■"-0.0 
49+1.3 
^•^-1.2 

o^o-o:^ 
11.0!" 

0.0+'-^ 



-0.6 
+12.0 



4.9! 
0.2-' 



-0.0 

0^0!^:^ 



3.5![- 
0.0+2: 



59.2 
13.1 



9.0 
0.0 



-3.7 
-1.1 
-0.0 



0.0 

10.2+3-8 

2 0+'°-'5 

^■"-0.0 

4.6+I-' 



25.5+ 



-0.0 

;+5.0 

-5.2 

6.8+8-^ 

12.4!4-4 

33.4+8:6 



19.1 

3.1 



10 



15.1 
26.8 



5.8- 



-4.2 



0.1 



-0.6 



+i¥o 

+i.-| 

10.9![4 

9 0+3.4 
^■'^-o.o 

7 0+2-2 
'■"-1.3 

12 5+5° 

j+2.2 
-0.0 

^■2l:^ 

+2.8 

-2.3 

+2.3 

3.4 

28.0!^ 

10.8!i.^2 

7 6+'-^ 
' ■"-o.o 

20.6!|-^ 

i6.i!2:4 

11 4+3-3 
^^■^-1.1 

12.5+" 
70 1+'^-'' 

'"■^-13.0 

7 9+'-^ 
'■^-2.3 

43.2!5;6 

1 3 7+6-'' 
^•'■'-0.0 

56 9+3-8 
-'"■^-11.5 

T T+4.9 

•'■•'-0.0 

8.3!4-6 

43 4+8'.2 

^■'■^-12.1 

7 4+10.1 

-2-.,3 , 



11.9 
0.0 



30.3 
3.8: 
13.3 
19.7 



3.2 
+9.9 

0.0 
+12.2 

15.5 

+2.5 

i 

-7.4 

9 1+3-3 
^■^-1.6 

0+5-6 

"■"-0.0 

3 2+'-2 
■'■^-0.0 

11 3+'-'' 

34 4+6.1 
■'^■^-4.9 
q •3+2.2 

9 0+'-5 
^■"-2.1 
44+11.7 
^■^-0.0 

1 6+3° 

^■"-0.0 

I T q+27.4 
^•'■^-O.Q 

2 5+1.6 
-1.6 

8.7!|:2 

6 1+" 

1 1+2-' 

0+6-' 

"■"-0.0 
j+1.9 

7+6:0 

'-2.2 
1+2.8 
'-0.0 

6.11:^ 

II 8+'°-' 

0+^'-'^ 

"■"-0.0 

0+4-3 
"■"-0.0 

15 7+6-2 

^-'■'-8.3 
+9.9 
12.9 



3.6 



1. 

9.7! 
2.3+ 



o 



M 
xj 



o 
3 
•a 
o 



n 

en 






24.5 



Table A.l. continued. 



to 

o 



0.1 /im 



l.SyLfm 



6.0;um 



Starname 



Am 



10.2!f« 
46.5+12-2 

54.9!^ 
21.2!5-5 

7 9+4.6 
'■^-0.0 

48 5+"-2 

2 0+"-^ 
21 4+2'' 
19.7!ji 
20 1+^-2 

18 4+7.3 

24.6!^:« 

25 6+5'^-^ 

1 5+2-^ 

23.6!^-4 

0+3-* 
q Q+5.2 
^■^-2.0 

12.1+31 
1 4+2 2 

6 1+2-2 

1(C Q + 11.5 

37-2!27 
24.5!t;S 
66.7!^2^o 
14 1+'-' 

i4.8!3-;5 

10.1-3 
4 6+2-" 

^■"-0.6 
o Q+130.1 
°-^-0.0 

7 7+5.9 

'■ -1.9 
o 7+12.2 

-%% 

+2.6 

-+°6^ 

;^:^ 

-3.0 



Ens 

__%_ 

0+2-5 
3 6+'-^ 

16.7+2-^ 

2.3+2-3 

6+^-'^ 

-Oft 
9 9+3.0 

1 2+''-^ 

^■^-0.0 
Q Q+6.4 
^■^-5.0 

0.0+°-^ 



For 

__%_ 
0+'-' 

"■"-0.0 
0+1-8 
"■"-0.0 

0+'-^ 

"■"-0.0 

-21.2 
-18.1 
1.6 
0.0 
1.5 



Sil 

% 

TTT" 
-0.0 

1 3+3-8 

^■■^-0.0 

0.0+1-1 



Am 



^3l^ 
11.4!|^ 

7 T+4.3 

'■-"-o.o 
25.8!t^ 

■^^■"-8.2 

63 9+5 5 
5 6+1-3 
69.9!^-4 
12.0!4-^ 
67.3+4-3 
55.1!5| 
44 3+^-^ 
28.0+4-^ 

68.6!!?i 
58.6+;^-^ 

+2.5 
-1.2 
141.2 
0.0 
17.6 



Ens 

% 
6 2+''-'* 

U./._j (J 

97 Q+16.9 
^'■^-13.6 

5.2+2-3 
8.0+5:: 



For 

% 

"■"-0.0 

0+2-'' 

"■"-0.0 

4.2+3-3 



Sil 

% 

oo+i-'* 

"■"-0.0 

0+2-** 

"■"-1.8 

0+1-° 

"■"-0.0 

0.0 
0.0_ 
0+*^ 

"■"-0.0 

0+1-1 

"■^-3.0 

0+1-5 

"■"-0.0 

1+1-*^ 

"■^-0.0 
"■"-0.0 

3.9!i-3 
0^0!g 

o^oil 

"■"-Q.Q 
"■"-0.0 

0+1-3 

"■"-0.0 
q q+320.1 
^■^-0.0 

0.7^5-0 

0+3:° 

"•"-0.0 

6+3-3 
"■"-0.0 

6+^-5 

"■"-0.0 

0.0+1-4 



Am 

% 



Ens For 



Sil 



'0 



SSTc2dJ162245.4-243124 

SSTc2dJ162715. 1-245139 

SSTc2dJ162816.7-240514 

SSTc2d J 182850.2+00950 

SSTc2d J 182900.9+02931 

SSTc2dJ182909.8+03446 

SXCha 

Szl02 

Sz50 

Sz52 

Sz62 

Sz73 

Sz76 

Sz96 

TWCha 

V710Tau 

VWCha 

VZCha 

WXCha 

XXCha 



1.9 



1.6 
19.6 
16.6 



0.0^ 

3.V 



n 



o^o!^:l! 
1 3+3-5 

^■-"-0.0 

1 3+'-' 
^■--"-o.o 

4 5+3-' 
0+5"*'-'^ 

"■"-0.0 

2^4!ii 
16.9li 

3 9+'-2 
-'■^-0.8 

0+2-3 
"■"-0.0 

2 4+1.2 

-'■-'-3.0 

7+1.1 
"■ ' -0.0 

T 7+3.6 
--'■ -0.0 

4 6+1-2 

^■"-2.4 

15.3!3-2 

9 7+4.2 

22+?:§ 
^■^-1.1 

T 7+2.4 
--'■ -0.9 

2.3+'-- 

10 5 +"^'^^-5 
^"■-"-0.0 

2 4+'-'^ 
-Q-8, 
11.3+ -,1 



1.7 
7.8 
2.4 



37.4: 

0.0; 
0.1 



-0.0 

-0.0 
0.0+0-5 



0.0+ 



0.0+ 
0.0; 
0.2^ 
0.0 



-0.0 
+3.9 
-6.5 



n 
^-0.0 

0+3-'' 
"■"-0.0 

"■"-0.0 

1 7+4-5 

^■'-1.2 
4C + 163.4 

"■"-0.0 
9 1 +8.0 
^■-^-3.2 

0+'-' 

"■"-0.0 

6 3+4-** 

H.3 



0.0 
0.0 



:?:8 

0.0 

+0.4 



3.9 

+2.2 
-1.4 



0^0!o.o 

o^o!ii 

2 1+'-** 
^■^-1.2 

3 i+i-'' 

0+4-1 

"■"-Q.O 
"■"-0.0 

0+2-6 

"■"-0.0 

0+2-' 

"■"-0.0 

3+°-4 

A ■7+46.6 
^■'-0.0 

2 6+3-3 
-9-1 

0+4-4 
"■"-0.0 

"■"-0.0 

0+3-^ 

"■"-0.0 

0+4-' 

"■"-0.0 



+5.7 
0.0 

+ 1.4 
1.3 

+ 1.0 
0.0 

+0.4 



+§:? 



6.3 
3.6 
0.0 
0.0 

5 1+^:^ 

-'•-'-2.7 
9 4+0.7 

^■^-0.7 
3 + 1.7 
"■•'-0.0 
9 0+0.6 

0.0_ 

"■"-0.0 

0+4-3 
"■"-0.0 

"■"-0.0 

9 T+4.6 
^■■-'-0.0 

1 7+89.9 

1 2+1-1 
-0.5 

0+6-° 

52;?:^ 

2+6-5 

"•^-0.3 

0+1-1 

j.z._22 

; + 1.0 

-0.0 

+8.2 

^-7.8 

0+1-1 

"■"-0.0 

"■"-0.0 

1 1+2-3 
^■^-0.0 

1 6+2-'' 

^■"-0.0 

0+"-3 

"■"-0,0, 

6 1+54-3 
"■^-0.0 

18.0!'" 

1 6+1-'' 

9 7+4.8 
^^ -0.6 



0.5^ 
14.3+ 



5.5 



-u 



52.8 

42 6+'4 

^^■"-14.6 

24 1+4-1 

^^■^-4.6 
CO c+22.0 
JO-J_14.8 

46.2!i'!,<' 
44.7!5-2 

27.9!^:3 
38 4+6-" 

-2.8 

68.0^4-^ 
1 7 9+4-2 

^'■^-3.8 

o^o!^:^ 

28.6!4-i 

23.8!23^-' 
+9.9 

iy 

+^6 



23.8+ 



55.0 
18.8 
56.1 



22 0+**:2 
^^■"-7.0 



-0.0 



o^8l:;! 

9 q+2.6 

^■^-2.2 

3 6+2-'* 
42+?:^ 
8.4+g 

+ 15.0 



28.3 
5.6 
4.1 
6.5 
6.0 



-18.1 
1.4 
1.7 
1.9 



i 

-2.3 

4 1+4.0 
^■^-0.0 

"■"-0.0 

io.o!2-3 

13 0+426-3 

^•'■"-0.0 

c q+4.8 

-"•^-4.4 

0+10.1 

"■"-0.0 

7 0+3.4 

'■■'-4.2 

10.8+^34 

3.2+2-1 
7.9+^:^ 



2.4 

+2.3 



7.9 

0^0_o.o 

5.6!|| 

8.4!3.| 
13.3^^42 

T 9 + l.'9 
•'■^-1.8 
0.0+''4-4 



+°4°9 
-^6.12 



0-0-0.0 

3^81:^ 

^+7.1 
-0.0 
>+9.6 
'-3.6 



0.0! 

14.2+ 



4.4; 
12.0 



2.2 

+1.8 

4.4 

0+2-3 

19+i§ 

3+2-2 
"■■'-0.0 

7+3.0 
•'■'-2.3 

1 7+2-5 
^■'-1.4 

0+1-0 

"■"-0.0 

0+"-** 

49+3:9 

-2-5 

0+"-' 

Ol+§:^ 
"■^-0.0 

3.5+!-5 



8.9 
0.0 



3.4 

+63.2 
0.0 
+3.6 



7.0_ 

0+^:^ 
"•"-0.0 

19 3+3 5 

-Ai 
0+" 

"•"-0.0 
4.3^3-1 

0+1:^ 

"•"-0.0 

0+1-5 
"•"-0.0 

1 5+1-8 

^•-"-0.0 

8.7!i-8 

2 3+3:3 
^••'-0.0 

8 1+3.5 

0+^:^ 
"•"-0.0 

18. 23-3 
22 0+23-5 

^^•"-7.0 

4 7+2.0 

^•'-1.6 

13 6+112-3 

8.0+^-^ 

13.4_ 

2 9+^-^ 

56+^:^ 

j.u 3Q 

115 + 11-6 
^^•-'-4.0 



AIa 



0.0 
+0.9 



7+3 



0.0! 
17.2+ 



4 7+ 

^- ' -2.5 

-1+1.6 

-0.0 

1+5.4 

--9.4 

0+°-^ 

10+?:8 
^•"-0.0 
4+3-^ 
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Table A.2. Dust temperatures, continuum parameters (blackbody temperature 
and offset Ov.i) derived using tiie "B2C" procedure for 60 objects from our sam- 
ple. 



Dust temperatures 



Continuum 



Goodness 



Starname 



(K) 



228.9+:;i' 
207.3+f^ 

238.0!^^ 

192.8+2-^ 
341.4+[0] 

252.6+^Y 
226.9!^| 

226.0+^1 

224. 1+'*-' 

231.8:f^ 

188.4+}] 

1 84.9+' 'I 

211.8+P 

225.9+1' 

194.1+!-" 



208.2 



Tcold 

(K) 



^ cont 

(K) 



58.2+o;g 
79 5+'-^ 
98.3+1^ 

Tc 9+4.9 

/J.Z,_j 

80.9!'| 
88.8+;;o 

102.9!^-^ 

91 0+^'^ 
aq 4+12.9 

87.1![^ 
85.2!;^ 

53 1+''' 
65 2+'-^ 
107.3!4-2 
83.4!2-6 
74 0+^-^ 

'^■^-2.0 

82.3+^-0 
48.9!i| 
55.6!'-^ 
92.3!| 

59 7+'-'' 

103.0!^-^ 

80.8!|-^ 

60.6!J-^ 

101.0!^^ 

13 1+'-^ 

73.5+2] 

91.4!« 

76.6!'| 

96 8+^-2 

61.6+'-^ 

60.2+2* 

58.6!2:5 

1+2:2 
■^-0.7 
90 3+2-^ 

101.5+I040 

59.0+'-5 



121.6"^ 

91A^_\\ 

107.6!gj 
+0.0 

:|:5 

-'-0.0 
85.2+":^ 

119.9+[-^ 
+1.3 

-6a2 
+0.8 

1.4 

94.0!ji 
83.5:«-^ 

81-895 
95.6!'f 

90.8+[-^ 

112.5+H 



(X 100) 



x; 
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TUT" 
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BD+31634 

C7-11 

CK4 
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HTLup 
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Harol-4 
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ISO-ChaII54 
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25.6 
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46.3 
10.6 
8.4 
39.6 
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32.6 
37.5 
31.0 
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34.3 
11.4 
25.3 
92.8 
14.5 
17.1 
18.4 
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29.1 
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11.0 
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26.7 
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7.6 
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Dust temperatures 


Continuum 


Goodness 


Starname 


^ warm 

(K) 


Tcold 

(K) 


^ cont 

(K) 


(X 100) 


^1 


Dof 


TWCha 

V710Tau 

VWCha 

VZCha 

WXCha 

XXCha 


260.3+«« 
234.3!« 
209.0+2-2 

233.6+^^5 
239.5+|4 

223 3+1'''* 

^^■'■•'-10.8 


91.2+'-'' 
19.0+!'? 

38.6+[^^ 
82.5+^-^ 
102.2+4-^ 


72 0+"** 
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i39.4i: 

100.4+1° 
qq q+o.O 


34 7+"-^' 
10.0+g 

"■^-0,2 

21 3+°-^ 

8 A+0-5 


10.8 
20.5 
12.2 
9.7 
11.6 
7.4 


218 
221 
232 
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244 
245 
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Fig. A.6. Fits to the 60 objects using the B2C procedure. 
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Fig. A.7. Continued 
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Fig. A.8. Continued 
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